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A microscopically motivated theory of glassy dynamics based on an underlying random first or- 
der transition is developed to explain the magnitude of free energy barriers for glassy relaxation. A 
O ' variety of empirical correlations embodied in the concept of liquid "fragility" are shown to be quan- 

titatively explained by such a model. The near universality of a Lindemann ratio characterizing the 
maximal amplitude of thermal vibrations within an amorphous minimum explains the variation of 
fragility with a liquid's configurational heat capacity density. Furthermore the numerical prefactor 
of this correlation is well approximated by the microscopic calculation. The size of heterogeneous 
reconfiguring regions in a viscous liquid is inferred and the correlation of nonexponentiality of relax- 
ation with fragility is qualitatively explained. Thus the wide variety of kinetic behavior in liquids 
of quite disparate chemical nature reflects quantitative rather than qualitative differences in their 
S . energy landscapes. 
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It is believed all classical fluids could form glasses if cooled sufficiently fast so as to avoid crystallization. Central 
to glass formation is a dramatic slowing of molecular motions on cooling the liquid. The existence of a universal 
description of glass transitions is suggested by empirical observations connecting deviations from the Arrhenius law 
| for the slowing of rates, nonexponential relaxations in the super cooled liquid state and the behavior of thermody- 
£^ ■ namic properties on cooling(l). Quantitative differences in behavior of different substances sometimes obscure this 
universality. This has led to a classification of liquids into "fragile" ones like o-terphenyl, having the most dramatic 
deviations from the Arrhenius law, and into "strong" ones like pure SiC"2 where the Arrhenius equation works well(l). 
In this paper, we show how the fragile versus strong behavior of liquids can be understood within a microscopically 
motivated theory based on the idea that glassy dynamics is caused by an underlying thermodynamic, ideal "random 
first order" transition(2-7). 

The notion that a random first order transition lies at the heart of glass formation received its early theoretical 
support from the remarkable confluence of approximate microscopic theories of the liquid glass transition(8-10) and 
■ the behavior of a large class of exactly solvable statistical mechanical models of spin glasses with quenched disor- 
der(ll). Two closely connected theories of the liquid glass transition suggest features similar to first order transitions. 
One of these, the so-called mode-mode coupling theory(8,12), focuses on the feedback between the slow fluctuations 
of fluid density in a molecule's environment on the motion of that molecule. This theory predicts a sharp transition 
in the dynamics as well as a characteristic behavior of the time correlation functions near the predicted transition. 
Mossbauer effect(13) and neutron scattering(14) are roughly consistent with these precursor phenomena. At temper- 
atures below the transition, mode coupling theory predicts the freezing of the liquid's configuration near to a given 
random configuration; i.e., there is broken ergodicity. Another approach to the glass transition directly addresses 
broken ergodicity by investigating the stability of a frozen density wave using either self-consistent phonon theory(9) 
or the density functional theory of liquids, applying them to aperiodic structures (10,15). The mode coupling, self- 
consistent phonon and density functional approaches all predict that there is a Lindemann criterion for the stability 
of an aperiodic density wave: just as for a periodic crystalline solid, thermal vibrations cannot yield a root mean 
square displacement of particles from their fiducial location exceeding roughly one tenth of the interparticle spacing. 
The precise value of the Lindemann ratio only weakly depends on the detailed intermolecular forces. The predicted 
Lindemann ratio corresponds well to the experimentally measured magnitude of the intermediate time plateau in the 
structure function measured by neutrons(14). A finite Lindemann ratio would be consistent with a first order phase 
transition, but glass transitions in the laboratory do not show a latent heat as ordinary first order transitions do. 
This lack of latent heat is explained by the existence of the large number of aperiodic structures that may be frozen 
in at a glass transition in contrast to the unique periodic structure formed in ordinary crystallization. Many exactly 
solvable models of disordered magnetic systems have been shown to exhibit freezing into many structures(ll,16,17). 
The major class of these also show a first order jump in a locally defined order parameter without any latent heat. 
This defines what has been called a "random first order" transition. Unlike Ising spin glasses, these models possess 
no symmetry between local states but have long range, quenched random interactions. Such systems include Potts 
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spin glasses(ll), p-spin glasses(17), and the elegantly solved Random Energy Model(16). There are further parallels 
between these systems and the phenomenology of glass forming liquids, most notably both glass forming liquids and 
these models exhibit a Kauzmann entropy crisis, i.e., the configurational entropy vanishes at a finite temperature 
above absolute zero(18). This crisis would define an underlying ideal glass transition. Whether the crisis for liquids 
would be avoided in some way at lower temperature than measurements have been made is controversial and is of 
limited relevance to describing the observed behavior using the analogy. In the exactly solvable statistical mechanical 
models, a dynamic transition occurs at a high temperature T^ coincident with mode coupling and stability analyses, 
but the thermodynamic transition does not occur until at a lower temperature, Tk the configurational entropy of 
different frozen solutions vanishes(4). The idea then is that the glassy dynamics in the measured temperature range is 
governed by the approach to an ideal glass transition described like that shown in the exactly solved models. There are 
two seeming differences between the exactly solved models and the situation for the liquid-glass transition. First, in 
liquids there is no quenched randomness; it must be self-generated. Second, while the models have infinite range forces, 
interactions in liquids are of finite range. The absence of quenched randomness has been addressed by exhibiting sev- 
eral mean field models without quenched randomness that do generate randomness internally(19-21). Also the formal 
statistical mechanical tools used for quenched random Hamiltonians, e.g., the replica technique, have been shown to 
be applicable to atomic fluid systems with self-generated randomncss(7). Furthermore, computer simulations of fluid 
glass transitions show replica symmetry breaking like a random first order Transition(22). The consequences of finite 
range interactions are more important. The finite range causes the dynamic transition at T^, like a spinodal of an 
ordinary first order transition to be smeared out. It becomes a crossover to activated dynamics. Below Ta, motions 
in the finite range system can still occur that involve the rearrangement of large regions of the liquid. The transition 
to such collective activated events in liquids has been confirmed in simulations(23,24). The events are driven by the 
configurational entropy. For finite range systems approaching a random first order transition, an "entropic droplet" 
scaling argument for the activation barriers naturally explains the non-Arrhenius transport behavior and leads to the 
Vogel-Fulcher law(4,5). The idea that configurational entropy is needed for motions in glasses predates the random 
first order transition theory and was described by Adam and Gibbs(25). The older argument is really quite different 
from the random first order transition theory, since it provides no explanation for how a rearranging unit's activation 
energy is related to the microscopic forces. Here we show how the near universality of the Lindemann ratio explains 
the connection between barrier heights and thermodynamics for liquids of varying fragility. 

The naive density functional approach used to obtain the Lindemann criterion for vitrification allows an estimate 
for the free energy of dynamic rearrangements. The density functional(10,26) assesses the cost of forming any density 
wave by breaking the free energy into an entropic localization penalty and an interaction term. 

F = J /(p(r))d 3 r = k B T J d 3 rp(r)[lnp(r) - 1] + J J d 3 rd 3 r'(p(r) - p )c(r - r')(p(r') - p ), (1) 

where po is the mean density. The localization cost is the same as for a perfect gas while the interaction term 
involves the direct correlation function of the liquid, a renormalized form of the bare interaction potential. The 
direct correlation function is determined by the condition that the functional gives small fluctuations in density 
reproducing the static liquid structure factor. Higher order terms in the density can also be included. In the 
frozen aperiodic state the density wave is decomposed into a sum of Gaussians centered around random lattice sites, 
p(r) = ^V(^) 3 / 2 exp(— a(r — r^) 2 ), where a represents the effective local spring constant that determines the rms 
displacement from the fiducial lattice site. The localization sites are {r^}. For large a, the densities around different 
sites overlap weakly giving 

| = k B T[ 3 - ln(^° ) - | + 1 ^ dW(p(r) Pa )c(r r')(p(r>) p Q ), (2) 

where N is the total number of particles and ro is the mean lattice spacing. We can take po^o = 1- For small a, F/N 
reduces to the perfect gas value. 

A similar free energy expression is obtained from self-contained phonon theory where the direct correlation function 
is replaced by the Mayer function / = e -/ 3 "( r ) — i for hard potentials(lO) or by the potential itself (7) . The free 
energy varies with the particular arrangement of sites {r^}, but assuming a is the constant, the mean free energy of 
aperiodic structures is plotted in Figure 1(a). The lowest value of a for which a secondary minimum occurs is given 
by the Lindemann value oil . This minimum representing the frozen wave is higher in free energy than the a = fluid 
phase. For the exactly solvable random first order transitions the excess free energy of the frozen solution is known to 
equal the configurational entropy of possible mean field solutions, TS C . For the fluid system in addition to the a = 
and a ~ oil uniform stationary solutions of the variational equation SF = 0, there are saddle points representing 
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droplet configurations in which a region of low a « forms in the midst of a given large a solution (See Figure 1(b)). 
This saddle point is a transition state for reconfiguring the frozen density wave. Within the melted region there is 
a multiplicity of states corresponding to other aperiodic arrangements of the atoms. Much below Ta, the interface 
should be quite sharp, that is, in a single atomic layer a changes from a value near to near zero. Close to the 
transition should be smoother with a slowly varying over many atomic layers. In both cases, there will arise a surface 
tension a reflecting the deviation of a in the layers with a different from the bulk free energy minima values. The 
density functional expression for the droplet free energy then is given as a function of the radius of the droplet much 
as in conventional nucleation, 

F{r) = -^nTs c r 3 + inar 2 . (3) 

Here s c is the configurational entropy density. 

The maximum of F(r) gives a reconfiguration barrier AF$ = ^-wa 3 /(Ts c ) 2 . A detailed calculation of this barrier 
for a specific glassy system, the random heteropolymer has been given by Takada and Wolyncs(27). This naive droplet 
result(4,28) differs from the Adam Gibbs suggestion AF$ = s*A/i/s c (25), where A/x is a bulk "activation energy" 
per particle and s* is the "critical configurational entropy" taken to be usually fcsln2. The AG formula is not the 
result of a self-contained microscopic calculation but assumes the free energy cost of dynamically reconfiguring a 
region is independently given from the free energies that determine the low energy structures themselves. There is 
no apparent reason to assume Afj, a constant for different substances. On the other hand, the modern random first 
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order transition theory does suggest universality for a based on the universality of the Lindemann ratio a L /tq. We 
see this in the following way: assuming a sharp interface between the localized and delocalized regions, the energy 
associated with the interface should be one-half of the interaction part of the free energy in the bulk stable phase. 
Therefore, a = ^ro[|nfes ln(^-) — s c (T)], where n is the density of particles. Since the localization part of the free 
energy depends only logarithmically on a, we can replace a by its minimum value az, which it achieves at Ta- Near 
Tk, on the other hand, we can neglect the configurational entropy part of the expression. For temperatures between 
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Ta and Tk, the errors of making these two approximations largely cancel. This gives a = ^titqUbT ln( ) = cfq 
as an approximation for temperatures much below T^. The universality of the Lindemann ratio a L /tq means 
a/nroksT is nearly universal and therefore that AF$ increases more rapidly with cooling for substances with a 
large configurational heat capacity. This explains the empirical correlation that strong liquids with nearly Arrhcnius 
rate slowing have small excess heat capacities contrasting with fragile liquids having large excess heat capacity with 
dramatically non-Arrhenius slowing. Near T^ the interface broadens and the sharp interface approximation breaks 
down. A gradient expansion of the free energy as a function of a yields a surface energy vanishing near T^. The 
universal value of <jq is only approximate. For given substance, the remaining temperature dependence of a from the 
broadening of the interface implies that the apparent fragility of liquids measured at high temperature should be larger 
than that measured at low temperature, as noted by Angell in his detailed survey of viscosity data(29). Similarly 
we note that a depends on the density and therefore the pressure. Thus although a kinetic glass transition defined 
by a specific numerical barrier height or fiducial relaxation time will be largely a function of the configurational 
entropy density there will be another explicit but weak thermodynamic dependence on pressure too. Consistent 
with Nicuwcnhuizcn's recent analysis of the dynamic effects on glass transitions caused by pressure and temperature 
change(30), this could explain the mild deviation of the Prigogine-deFay ratio from 1. 

While the simple density functional calculation explains qualitatively the fragility/heat capacity density correlation, 
viscosity data are more consistent with an s" 1 scaling for the free energy of activation (like that suggested by Adam 
and Gibbs(25) ) rather than the s~ 2 behavior predicted from the simple density functional theory. The scaling 
theory of the cntropic droplet formulation already accounts for this observation(5). The modification comes from 
the complexity of the interface between aperiodic crystalline minima(5). Correct scaling near Tk is restored by the 
wetting of droplets corresponding with one particular density wave, by a surface coating corresponding to a different 
aperiodic arrangement. This acts to lower the surface energy much like what happens in the random field Ising 
model(31). Wetting for a random system gives a surface tension that depends on the radius of the drop. This r 
dependent energy yields sj 1 scaling when the thermodynamic critical exponents for the random first order transition 
are used. We now reprise this argument based on a similar one for the random field Ising magnet (31) in Figure 1 (b). 

The wetting argument leads to a differential renormalization group equation for cr(r), 

a^da = -(4-V3 - 4- 4 / 3 )(TV^A^) 4/3 r- 5 / 3 dr, (4) 

where Ac p is the heat capacity jump per unit volume. This renormalization group equation is integrated outward 
from r where the short range value is set by the naive density functional theory without wetting discussed earlier, 
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(To- Between Tk and Ta, a(r) vanishes at large distance and is only finite below Tk- Using this boundary condition, 
the solution of the renormalization equation for cr(r) at Tk is then 

a(r)=a Cf)^. (5) 

When this is substituted into the expression for F(r), one finds that the maximum gives a barrier, AF$ which now 
varies inversely to the first power of the configurational entropy density; i.e., the Vogel-Fulcher scaling. We find a 
simple expression for the activation barrier: 

Ts c TAc p T — Tk T — Tk V 1 

The coefficient D in this expression has been called the liquid's fragility, which has the expression 

D= 2 l„VpL ln >^A. (7) 
16 Ac p ire y ' 

— 1/2 

Based on the Lindemann ratio universality, the root mean square displacement, a L is taken as 0.1r , the hard 
sphere value, so that D can be expressed in terms of the heat capacity jump per mole, Ac p , 

D = 32i?/Ac p , (8) 

where R = 8.31 J mole -1 K _1 . The value of D depends on the heat capacity jump per mole which varies greatly 
from substance to substance and is far from being universal. In Figure 2 we plot the D predicted from this theory 
versus the inverse of the configurational heat capacity for several glass forming liquids. The straight line is given by 
Equation (8). Superimposed on the plot are the experimental values of the D. The agreement is excellent. We see 
that the magnitude of the activation barriers for rearrangement of the viscous liquid depends on the difference in 
temperature from Tk, on universal microscopic parameters connected with the Lindemann ratio and on the excess 
heat capacity connected with configurational excitations. 

A hallmark of the random first order transition theory of glass dynamics is the dynamic heterogeneity required 
to explain the growing barriers upon cooling. After combining Equation 3 and 5 with our expression of ao, a little 
algebra shows the characteristic size of a rearranging region is ^ = 2( 2 air i ) 2 ^{ t-t k Y^- This can be expressed 
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Vogel-Fulcher law. This is plotted in Figure 3. 

The kinetic laboratory glass transition occurs when molecular slowing gives relaxations in the hours range, i.e., 
^- = 10 17 . Thus at T g , (^) 4.5, a rather modest size. We also note the universality of cr suggests that s c is nearly 
the same for all substances at the laboratory glass transition. Roughly 90 molecules are involved in a rearranging 
unit according to the random first order transition theory at the conventionally defined glass transition temperature. 
The rearranging unit according to the Adam-Gibbs argument is a region just capable of having two states, therefore 
i^ 2 ') = ( fl '° 2 ) 1 ^- £ag grows slowly as Tk is approached in contrast to the random first order transition theory. 
The Adam-Gibbs argument gives rearranging units with at most 10 molecules near T g for the most fragile substances. 
It is clearly very ambiguous to have such small "cooperative" units. Recent observations of structural heterogeneities 
are inconsistent with units of the small size predicted by Adam-Gibbs, but arc in harmony with the estimates of 
the random first order transition entropic droplet picture(32). A single size does not characterize the viscous liquid 
completely. The random first order transition entropic droplet picture actually leads to a "mosaic" structure of the 
supercooled liquid(33) with cooperative regions only somewhat larger than the critical droplet size £. These regions 
fluctuate in size and therefore have different flipping rates because of the configurational entropy fluctuations whose 
magnitude depends on the configurational heat capacity density jump Ac p and the volume of the rearranging region, 
AS C = ^/fcsAcp^ 3 . At T g both strong and fragile liquids have roughly the same absolute scale for their mosaic 
structures, i.e., £ is nearly universal at the laboratory glass temperature. It follows that the range of activation 
barriers is smaller for strong than for fragile liquids because of their smaller Ac p . This is in accord with the observed 
correlation between growing nonexponentiality of relaxation with growing fragility (34). 

We conclude that the wide variety of kinetic behavior seen in liquids reflects quantitative rather than qualitative 
differences in their energy landscape. Furthermore the random first order transition approach coupled with microscopic 
considerations about the stability of aperiodic structures can account semi-quantitatively for these differences. 
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Figure 1 a, Free energy as a function of order parameter a. Right below Ta, a second minimum emerges around 
a w ax, which corresponds to a glassy state. The free energy difference between the liquid and glass state is TS C (T), 
which approaches zero at the Kauzmann temperature Tk- b, An illustration of a liquid-like (multiconfiguration) 
droplet inside a glassy region corresponding to a single mean field minimum free energy configuration. The interface 
is wetted by suitable configurations to lower the surface energy. One considers an inhomogeneous situation with 
single minimum given by the density functional theory abutting another minimum as in a naive droplet solution 
with a radius of curvature r. Upon this surface, one erects a smaller droplet of one of the other solutions of the 
density functional theory as shown in the figure. The free energy of interpolating this wetting phase is given by 
SF = <7(r)r d_1 (-) 2 — 5s c r d / 2 (^) 1 / 2 . This additional free energy cost depends on the surface tension at the scale 
r, o~(r) and on the fluctuations in driving force for forming this smaller wetting droplet. In the Ising model, the 
fluctuations in driving force for these droplets arise from the random part of the magnetic field. For a random first 
order transition, the field fluctuation's role in the disordered magnet is played by the fluctuations of configurational 
entropy density. The magnitude of these fluctuations should be given by the usual Landau expression AS 2 = ksAC p 
where AC P is the configurational heat capacity of a region. The contribution to the free energy from the interpolating 
wetting droplet yields a change with size of the surface tension at size r, da (Equation 4). 



Figure 2 The fragility parameter I? as a function of the inverse heat capacity jump per mole. The glass formers 
chosen are those shown in Angell's review article (a. Angell, C. A. Formation of glasses from liquids and biopolymers. 
Science 267, f924-f935 (1995)). The solid line in the graph is calculated with random first order transition model 
(Equation (7) and (8)) based on the universality of the Lindcmann ratio and the points are from experiments. Data 
for fragility parameter D are found in a, b. Korus, J., Hempel, E., Beiner, M., Kahle, S. & Donth, E. Temperature 
dependence of a glass transition cooperativity. Acta Polymer 48, 369-378 (1997), c. Richert, R. & Angell, C. A. 
Dynamics of glass forming liquids V. On the link between molecular dynamics and configurational entropy. J. Chem. 
Phys. 109, 9016-9026 (1998) and references therein; for specific heat jump are found in d. Angell, C. A. & Smith, 
D. L. Test of the entropy basis of the Vogel-Tammann-Fulchcr equation: Dielectric relaxation of polyalcohols near T g . 
J. Phys. Chem. 86, 3845-3852 (1982), e. Angell, C. A. & Torcll, L. M. Short time structural relaxation process in 
liquids: Comparison of experimental and computer simulation glass transitions on picosecond time scales. J. Chem. 
Phys. 78, 937-945 (1983), and f. Torcll, L. M., Ziegler, D. C. & Angell, C. A. Short time relaxation processes in 
liquids from viscosity and light scattering studies in molten KC1 • 2BiCl3. J. Chem. Phys. 81, 5053-5058 (1984). The 
heat capacity jump is given as per mole "beads" (d and g. Wunderlich, B. Study of the change in specific heat of 
monomeric and polymeric glasses during the glass transition. J. Phys. Chem. 64, 1052-1056 (1960) or "mobile units" 
(h. Schulz, M. Energy landscape, minimum points, and non-Arrhenius behavior of supercooled liquids. Phys. Rev. B 
57, 11319-11333 (1998)). Generally speaking, "beads" are "rearrangeble elements in a relaxing liquid" (d). We have 
used the bead count of previous workers (g). The number of beads for Ge02 and ZnCkj is 3 (d), 6 for glycerol (d), 10 
for KC1- 2BiCl 3 (f), 12 for 3KN0 3 •2Ca(N0 3 ) 2 (e), 2 for m-fhiorotoluene (g and I. Chang, S. S. & Bestul, A. B. Heat 
capacity and thermodynamic properties of o-terphenyl crytal, glass, and liquid. J. Chem. Phys. 56, 503-516 (1972)), 
and 3 for o-terphenyl (i) . The "bead" count is a crude way of accounting for internal flexibility of the molecules (since 
the free energy functional is essentially that for a monatomic fluid) . To illustrate the robustness of the correlation we 
indicate also the values of o-terphenyl as shown by a star (*) if its internal flexibility is ignored and it is assigned a 
bead count of one. 



Figure 3 The correlation length £ (in the unit of lattice spacing tq) is shown as a function of relaxation time. 
The solid line is that predicted by random first order transition theory while the dashed line is the result of Adam- 
Gibbs theory(25) assuming Ac p = 51.8 J mole" 1 K _1 , the value for PVAC (polyvinyl acetate)(25). The Adam-Gibbs 
result weakly depends on fragility. The point (and its error bar) gives the only "direct" measurement by Spiess et 
al. on PVAC (32). Results of many "indirect" measurements on different glass formers around the glass transition 
temperature also fall in the range of the Spiess experiment (32). 
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Figure 1 a 




7 



Figure 1 b 
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Figure 2 
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